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Abstract 

Magnetic moments of baryons composed of light and strange quarks 
are computed for the first time through the only parameter of the 
model- string tension a. For example, fip = mpjc^^ ^i^- = — + 

— ^-^^)~^ ^ where rup is the proton mass, is the strange 
quark current mass, c = 0.957 - a constant which is caluculated in the 
paper. Resulting theoretical values differ from the experimental ones 
typically by about 10%. 

1 Introduction 

Recently in a stimulating analysis of the hyperon static properties H.J. Lip- 
kin displayed remarkably successful relations connecting strange and non- 
strange baryon magnetic moments and the corresponding quark masses Vu 
and Vg (in order to avoid similar notations for masses and magnetic moments 
we denote current quark masses by and dynamical, or constituent, quark 
masses by Vk)- To understand why these relations so well agree with the 
experiment and to get insight into the problems encountered in semileptonic 
decays of baryons one needs a dynamical approach which would enable 
one to express constituent masses of quarks and baryon magnetic moments 
in terms of a single QCD scale parameter. It is a purpose of the present letter 
to express the magnetic moments of baryons and constituent masses of the 
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corresponding quarks in terms of only one parameter - the string tension, 
and to demonstrate that our results are in line with the relations of Lipkin. 

Theoretical investigation of baryon magnetic moments (BMM) has a long 
history 0. In the constituent quark model (CQM) BMM are expressed 
through the values of constituent quark masses, which are input parameters 
(see also |Q for discussion). Among other approaches to the problem 
mention should be made of different versions of the bag model lattice 
calculations and the QCD sum rules 0. In the latter the BMM are con- 
nected to the values of chiral and gluonic condensates and to the quartic 
quark correlator. Although a lot of efforts has been undertaken along differ- 
ent lines, the theoretical predictions still differ from the experimental values 
(by 10-15% in the worst case ||). 

In all models however theoretical predictions are somewhat biased by 
the introduction of supplementary parameters in addition to the only one 
pertinent to QCD - the overall scale of the theory, which should be specified 
to make the QCD complete. In the final, ideal case this role is played by 
^QCD'i in our treatment, as well as in lattice QCD calculations, we take as 
this universal parameters the QCD string tension a, fixed in nature by the 
meson and baryon Regge slopes. 

The purpose of our letter is to calculate BMM through this single param- 
eter, 0", in the simplest possible approximation within the nonperturbative 



QCD approach, developed in P]-|14 



2 Relativistic 3g Green's function and effec- 
tive Hamiltonian 

The starting point of the approach is the Feynman-Schwinger (world-line) 
representation of the 3g Green's function 0, where the role of "time" pa- 
rameter along the path z^^\si) of z-th quark is the Fock-Schwinger proper 
time Sj, i = 1,2,3. One has P,pn]] 

G'(3.)(x,F) = J l[ds,Dz^;^e-''{Ws{X,Y)) (1) 

i=l 

where X;Y = x^^), x^^); yW, y(2), y(3)^ 

+ (^) *<). (2) 
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Here rrii is the current quark mass and the three-lobes Wilson loop is a 
product of three parallel transporters 



{Ws{X,Y)) = (n<i(^^^rO)ea,a,a3e.,M3- (3) 
1=1 

The standard approximation in the QCD string approach is the minimal area 
law for d^), which will be used in what follows 

{Ws) = exp{-aj2S,) (4) 

i=l 

where 5*^ is the minimal area of one loop. 

The next step is basic for our approach, and it allows finally to calculate 
the quark constituent masses z/j in terms of the quark current masses rrii, 
defined at the scale of IGeV. 

In this step one connects proper and real times (in the baryon cm. sys- 
tem) 

dsi = —T- (5) 

where t = z^\si),0 < t < T, is a common cm. time on the hypersurface 
t = const. The new entity, z/j(t) as will be seen, plays the role of the quark 
constituent mass and will be calculated through (j(and ctg when perturbative 
exchanges are taken into account). 

Considering the exponent in (^) as an action one can define the Hamilto- 



nian and go over to the representation [0,|]13 

G'(3'^) = J llDu,{t)Dh^'\t)e-^ (6) 

with 



i=l 



The final step in the approach ||lT|, |12| is the derivation of the cm. 
Hamiltonian containing Vi as parameters to be found from the condition of 
the Hamiltonian minimum. It has the following form 
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Here ^, r] are Jacobi coordinates defined as in [0, and r^*^) is the distance 
from tlie k-th quark to tlie string-junction position which we take below for 
simphcity coinciding with the c.m.point. In addition contains an arbitrary 
mass parameter m introduced to ensure correct dimensions, this parameter 
drops out from final expressions. Leaving technical details to the Appendixes, 
we now treat the Hamiltonian (P) using the hyperspherical formalism [O . 



3 Evaluation of quarks constituent masses 

Considering three quarks with equal masses and introducing the hyperradius 
p2 — ^2 _|_ Qj^g j^g^g jj-^ approximation of the lowest hyperspherical 



harmonic (which is known |jT6[ to yield accuracy of the eigenvalue around 
one percent) 



dp- 
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2K£„-»'(p)}x(p)=0, (9) 




The baryon mass M„(z/) is equal to (for equal quark masses) 

M^{u) = ^ + ^u + EM- (11) 

The crucial point is now the calculation of u, which is to be found from the 
minimum of Mnii^), as it is prescribed in the QCD string approach [P[]-||14[. 

At this point it is important to stress that we have changed from i/(t) 
depending on t on the trajectory in the path integral (|^) to the operator u to 



be found from momenta and coordinates in (^ as in and finally to the 
constant u to be found from the minimum of the mass M, as it was suggested 
in |T0[,[|11[], [|I^. The accuracy of this replacement was tested recently in 
to be around 5% or better for lowest levels. 

The equations defining the stationary points of M„ as function of u for 
equal masses is 

-Q^ 1-^(0. = (12) 

The generalization for baryon made of three quarks with different masses 
is straightforward. The perturbative gluon exchanges and spin-dependent 
terms can be selfconsistently included in the above picture W%- Including 
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the Coulomb term and passing to dimensionless quantities x, e„ and A defined 

as 

x = {2ub)^p, e^ = j^,. ^ = "^3i^^j ' ^^^^ 

where as is the strong coupling constant, one arrives at the following reduced 
equation 

x^-^---eri\)\x{x) = ^. (14) 
It is now a simple task to find eigenvalues en(A) of (0) either numerically, 




or analytically (see below). Then ([T^) would yield the following equation 
defining the quark dynamical mass v 




9 /757r2\'/' /m2 



16 \ 2 / \v 



— -1 =0. (15) 



d\ 

It turns out that numerical solution of (|T^) may be reproduced ana- 
lytically with the accuracy of (1-2)% provided one replaces the potential 
W{x) = X + d/x"^ ~ X/x in (p^ by oscillator potential near the stationary 
point W'{xo) = - see the Appendix A. Equation ([T5|) applied to the nu- 
cleon {m = 0) yields the dynamical mass of the light quark, and applied 
to Q~{m = irts) gives the strange quark mass v^. Before presenting these 
solutions we remind about spin-spin forces responsible e.g. for iV — A split- 
ting. Contrary to what might be naively expected the inclusion of spin-spin 
interaction considerably simplifies the problem due to remarkable cancella- 
tion of Coulomb and spin-spin contributions into the dynamical (constituent) 
quark mass - see the Appendix A. Therefore these terms should be kept only 
if one wishes to calculate the BMM with the accuracy much higher than 
10% in which case one should also take into account pion corrections, higher 
hyperspherical harmonics, etc. which is out of the scope of the present paper. 

Thus in order to determine the quark masses and eventually the baryon 
magnetic moments one needs only two parameters: the string tension a and 
the strange quark current mass (m, d current quark masses are set to zero). 
Present calculations were performed for 

a = O.lSGel/^ m, = 0.245GeV. (16) 

The string tension value (16) which is smaller than in the meson case is in line 



with baryon calculations by Capstick and Isgur [18|. A similar smaller value 
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of a is implied by recent lattice calculations by Bali [|T^. Since the value of the 
above parameters are allowed to vary within certain limits ||14||,||19||, [pO[ 



one 



can in principle formulate the inverse problem, namely express the BMM in 
line with the present work and then fit their experimental values to determine 
the optimal choice of a and rris. 

Consider first the case of a nucleon made of three quarks with zero current 
masses and equal dynamical masses u^. Keeping in mind cancellation of the 
Coulomb and spin-spin terms and thus setting in ( p!4D and (|l^) A = and 
making use of the oscillator approximation described in the Appendix A, one 
finds from ( jlS]) 




3 ■ 51/3 



1 + 



3/4 



= c 



^/a ~ 0.957v^ = 0.37GeV. (17) 



This result agrees with the exact solution of (|l^) at A = with the accuracy 
better than 1%. Similar procedure applied to Q~ baryon yields the strange 
quark dynamical mass z/g. From (p!5[) and ([T7| ) one gets 




15 ml 
32^ 



OAGGeV 



(18) 



for rris = 0.24:5GeV, and where the constant c is defined in ([T7|). Now we 
turn to baryon magnetic moments. 



4 Baryon magnetic moments 

The form (|I]) for G^^''^ does not take into account spins of quarks. When those 
are inserted, a new additive term appears in the exponent of (^), proportional 
to the external magnetic field B, namely A acquires the following term ||21| 
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6A = J: ' drkCCT^'^B = Y. -^^^dt, (19) 

where is the electric charge of the quark, (t^^'> is the corresponding spin 
operator, and the definition (^ of the constituent mass was used. 
Introducing the z-component of the magnetic moment operator 
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one can write the BMM as matrix elements 



(21) 



where is the eigenfunction of (^), and Uk is taken at the stationary point, 
given by (|l|). 

For the baryon wave function we shall take here the simplest approxima- 
tion, namely 

= ^«s'™™(r)V''^"'™(a, f)ij'' (color), (22) 

where ip{(T, /) is the spin- flavour part of the wave function. The form (0) ne- 
glects the nonsymmetric components in the coordinate iplr) and spin-flavour 
parts of wave function, which appear in the higher approximation of the hy- 
perspherical formalism and for lowest states contribute only few percent 
to the normalization [|15|. The spin- flavour functions for different baryons 
were known for a long time and are briefly outlined in the Appendix B. 
Using these functions it is a simple task to calculate the matrix element (pT]), 
e.g. the proton and neutron magnetic moments are given by 



rrir. 



Hp 



-1 -3/4 



3 ■ 51/3 



1 + 



3^5, 



~ 2.54/iAr, = -g/^p - -1-69/i 

(23) 

Magnetic moments of other baryons as well as new relations between them 
are obtained from (|1^), (|l3) and fH! 
In particular one has 



N- 



f^n- - -/ip 1 + - — 
4 c'a 
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(24) 



32cV2^ 

/is+(4cV + ml) = -2/2=0 (3cV + 2ml), (25) 

where terms of the order ml/c^a'^ were omitted in deriving the last relation. 

Results on the BMM are summarized in Table 1. As these results differ 
from the experimental values typically by only about 10% and are subjected 
to plentiful corrections (meson exchanges, higher harmonics, etc.), one may 
conclude that the outlined QCD approach is successful even in its simplest 
form. 

It is important to realize that the QCD string model used above is a fully 
relativistic string model for light current masses, and the "nonrelativistic" 
appearence of the Hamiltonian (|1|) is a consequence of the rigorous eibein 



formalism |23] which was introduced in the most general form in |24 
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The approach enables to investigate other electromagnetic properties of 



baryons: transition magnetic moments, polarizabilities [|26|, etc. 
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Table 1. Magnetic moments of baryons (in nuclear magnetons) com- 
puted using Eqs. (p!2D , (pOD , (pi]) in comparison with experimental data from 
PDG @ 



Baryon 


P 


n 


A 






s+ 




^0 


n- 


Present work 


2.54 


-1.69 


-0.69 


-0.90 


0.80 


2.48 


-0.63 


-1.49 


-2.04 


Experiment 


2.79 


-1.91 


-0.61 


-1.16 




2.46 


-0.65 


-1.25 


-2.02 
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Appendix A 
Eigenvalue equation for baryons 



The eigenvalue equation for baryons in hyperspherical basis [|T5| in its 
standard form is given by (Wj and reduced form with Coulomb term included 



by (0). Though (|T^ can be easily solved numerically, it is instructive to 
present the analytical solution resorting to the oscillator approximation near 
the stationary point W'{xq) = 0. This approximation is known to yield the 
accuracy of 1-2% [0. It will be demonstrated below that the Coulomb term 
which tends to increase the quark dynamical mass almost exactly cancel with 
the spin-spin interaction term. Therefore one can consider ( pTD at A = 0. 
Then the ground state energy is equal to W{xo) plus the first quantum 
correction uj/2, where a;^ = 2W"{xo). 
This yields 

Substitution of this result into (^) and (PTD-(P^) leads to the expression 
([T7| ) for the light quark current mass Uu- 

Next we demonstrate the cancellation of the Coulomb and spin-spin con- 
tributions into the quark mass. Again this conclusion results directly from 
numerical calculations but it is always preferable to present transparent es- 
timates. The derivative dEn/dv (see (|11|), (0), (0)) may be written as 



dEn 
dv 



e{\) 



d (2z/6)2/3 (2z/6)2/3 d\ 



du 2u 



2u du 



de 



dX 



(A.2) 



where the (-) sign stems from the fact that de/dX < 0. Expanding e{X) 
in Taylor series in A and keeping only linear term (the small parameter is 
A/e(0) ~ 1/4) one gets 

(2z/6)2/3 



dEn 
du 



-£(0) 1 + 



A 



de 



dX 



(A.3) 



The value of e{0) is given by ([A.l|) , the estimate of de/dX in the small A 
regime is strightforward, then recalling that according to ([TsD A oc (z/^/a)^/^ 
and solving the simple equations one obtains that due to Coulomb interaction 
Uu increases by ^ O.OSGeV^. This is confirmed by numerical solution of (|14|). 
The spin-spin interaction in baryon results in the shift of En equal to 



SE„ 



S7, 



16 ^ 

Cts / 

9 UiUj 



Sir, 



(A.4) 
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For proton the summation over {i,j) yields a factor —4/3, all three delta- 
functions smeared over infinitesimal regions [|16] are equal to each other and 



scale with u as z/^/^ ■ 5, where the constant 5 for nucleon has been with high 
accuracy computed by Green's function Monte Carlo method in [0. As a 
result spin-spin interaction leads to a contribution into dEn/dv proportional 
to z/~^/^. This in turn results in the decrease of the quark dynamical mass 
Vu by 0.035 GeV, i.e. the contributions from the Coulomb and spin-spin 
interaction into the quark mass almost exactly cancel each other. Thus we 
are led to a value Vu = 0.37GeV given by (p!7|). 

Appendix B 

Spin-flavour wave functions and BMM in impulse Approximation. 

As stated in the main text we have restricted the basis by considering 
only totally symmetric component of the baryon wave function in the coor- 
dinate space (see for the discussion of corrections to this approximation). 
Therefore the spin-flavour part of the wave function has to be symmetric 
too. Then the calculation of the BMM in "impulse" (additive) approxi- 
mation proceeds along the well trotted path [Q]. For example, the nucleon 
spin-flavour symmetric wave function entering into (^2]) is a combination 
^symmf^^^ /) = ((f" {a)x" (f) + ^' io-)x' (f)) / V^, whcrc prime and double prime 
denote mixed symmetry functions symmetric and antisymmetric with respect 
to the 1 — i> 2 permutation. The E~ spin-flavor wave function is obtained from 
that of neutron one by substitution of the w-quark by s-one and so on for 



other baryons. Due to the antisymmetry of the complete wave function (p2D 
the calculation of the matrix element (|21|) reduces to the averaging of the 
operator {S^u/ J^u + ^ss/ ^s) or {63d/ + ^ss/ J^s)- In this way one arrives at the 
well-known relations 

/in _ _2 flA _ /iE+ _ 8 _^ ^ - _1 + 

fip 3 ' /ip ?>i's ' /ip 9 9z/s ' /ip 9 9z/s ' 

^ - ^ _j_ J:^^ - i _ 1^ - _^ _ 1^ 

lip 9 9z/s ' /ip 9 9 fip 9 9 Vs 

Various corrections to impulse approximation have been discussed in the 
literature 01,11. 
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